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ABSTRACT 

We describe the characteristics and evolution of the magnetic field and chromospheric emission in an active region 
core observed by the Solar Optical Telescope (SOT) on Hinode. Consistent with previous studies, we find that the 
moss is unipolar, the spatial distribution of magnetic fiux evolves slowly, and that the magnetic field is only moderately 
inclined. We also show that the field line inclination and horizontal component are coherent, and that the magnetic 
field is mostly sheared in the inter-moss regions where the highest magnetic flux variability is seen. Using extrapolations 
from Spectropolarimeter (SP) magnetograms we show that the magnetic connectivity in the moss is different than in the 
quiet Sun because most of the magnetic field extends to significant coronal heights. The magnetic flux, field vector, and 
chromospheric emission in the moss also appear highly dynamic, but actually show only small scale variations in magnitude 
on time-scales longer than the cooling times for hydrodynamic loops computed from our extrapolations, suggesting high- 
frequency (continuous) heating events. Some evidence is found for flux (Call intensity) changes on the order of 100-200 
G (DN) on time-scales of 20-30 mins that could be taken as indicative of low-frequency heating. We flnd, however, that 
only a small fraction (10%) of our simulated loops would be expected to cool on these time-scales, and we do not find clear 
evidence that the flux changes consistently produce intensity changes in the chromosphere. Using observations from the 
EUV Imaging Spectrometer (EIS) we also determine that the filling factor in the moss is ~ 16%, consistent with previous 
studies and larger than the size of an SOT pixel. The magnetic flux and chromospheric intensity in most individual SOT 
pixels in the moss vary by less than ~ 20% and ~10%, respectively, on loop cooling time-scales. In view of the high 
energy requirements of the chromosphere, we suggest that these variations could be sufficient for the heating of 'warm' 
EUV loops, but that the high basal levels may be more important for powering the hot core loops rooted in the moss. The 
magnetic fleld and chromospheric emission appear to evolve gradually on spatial scales comparable to the cross-field scale 
of the fundamental coronal structures inferred from EIS measurements. 
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\^ Significant progress in solving the decades old coronal heating 
Qjfoblem could be made if one knew the duration and frequency 
Cbf heating events. Analysis of soft X-ray loops in the Yohkoh era 
^""Suggested that high temperature (3-5MK) coronal plasma could 
^e heated steadily (Porter & Klimchuk 1995; Kano & Tsuneta 
•"^MB). These loops are rooted in the "moss," which is found 
active region cores (Martens et al. 2000), and it has been 
^rgued that the lack of intensity variations there is indicative of 
■ - steady heating (Antiochos et al. 2003). Furthermore, hydrostatic 
modeling of whole active regions has been quite successful at 
reproducing the core emission from short hot loops (Schrijver 
et al. 2004; Warren & Winebarger 2006; Lundquist et al. 2008). 

Hydrostatic modeling, however, has greater difficulty repro- 
ducing the emission at lower temperatures. Warm (IMK) 
EUV loops observed by the Solar and Heliospheric Observatory 
(ySOHO, Domingo et al. 1995) and the Transition Region and 
Coronal Explorer {TRACE, Handy et al. 1999) have been found 
to be overdense compared to static equilibrium theory and per- 
sist far longer than expected loop cooling times (Lenz et al. 1999; 
Aschwanden et al. 2001; Winebarger et al. 2003). These obser- 
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vational features can be explained if coronal loops are bundles 
of unresolved threads that are heated impulsively (Aschwanden 
et al. 2000; Warren et al. 2002; Winebarger et al. 2003). 

Of course, most of the proposed coronal heating mechanisms 
are impulsive in nature (Klimchuk 2006), and the term 'steady 
heating' is usually taken to mean that the repetition time be- 
tween impulsive heating events is shorter than the time it takes 
for the loop to cool by conduction and then radiation. Loops 
are thus maintained at high temperatures and the emission is 
apparently steady. Many active regions appear to evolve slowly 
and loops are not often seen cooling in the core around the moss 
in these regions (Antiochos et al. 2003; Patsourakos & Klimchuk 
2008). In other cases, loops are clearly seen evolving and cool- 
ing (Ugarte-Urra et al. 2009), and it is not clear which type of 
heating is dominant. Still, it would be surprising if the structure 
and heating characteristics of warm loops and hot loops were 
fundamentally different. A challenge to current loop modeling 
is to draw these pictures together and understand how they can 
exhibit these apparently contradictory properties. 

The instruments on board the Htnode satellite (Kosugi et al. 
2007) are providing unprecedented observations of active regions 
in terms of temperature coverage, spatial, temporal, and spectral 
resolution. As such they are allowing us to probe the properties 
of moss at the bases of high temperature loops in active region 
cores in new detail. In two previous papers (Brooks & Warren 
2009; Warren et al. 2010, hereafter Paper 1 and Paper 2), we 
analyzed EUV Imaging Spectrometer (EIS) and X-ray Telescope 
(XRT) observations of an active region observed by Hinode in 
June 2007 (NOAA 10960) in order to study the time-scale of 
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energy release. Several lines of evidence led to the conclusion 
that the heating in the core of this region was effectively steady. 
First, soft X-ray and Fexii 195.119 A intensities in the moss 
were shown to vary by less than 15% over many hours. Second, 
from the Fexii 195.119 A line profiles, we measured Doppler 
and non-thermal velocities in the moss of ~ 3 km s~^ and ~ 26 
km s~^ on average, respectively. Taking into consideration the 
uncertainties, and also the fact that FeXII lines have been shown 
to be bluc-sliiftcd \w a few km s^^ in the quiet Sun (Peter & 
Judge 1999), the moss velocity measurement is consistent with 
zero. In addition, based on a comparison with several quiet Sun 
synoptic datascts, the measurement of the non-thermal velocity 
in the moss was shown to be no larger than the typical quiet 
Sun value of 25 km s^^. More importantly, neither quantity 
varied by more than 15% over many hours. Third, no evidence 
was found for co-spatial warm and hot emission, as would be 
expected from impulsive heating models that assume coronal 
loops have time to cool substantially between events. Finally, 
the observed moss intensities could be brought into agreement 
with hydrostatic simulations provided that local expansion of 
the magnetic field at the base of the corona was included in the 
numerical model. 

These results rule out the possibility of low frequency impul- 
sive heating of monolithic loops on the spatial scales resolved 
in the analysis. An alternative possibility, however, is that the 
heating is impulsive at low frequency on sub-resolution threads, 
as is thought to be the case in the warm overdense loops (Jim 
Klimchuk 2009, private communication). This picture may ap- 
pear like high frequency heating in the observational analysis of 
Paper 1. 

It is clear that we can already resolve the expected loop cooling 
times temporally, and even other transient activity on localized 
scales, e.g. short duration "Type II" spicules (de Pontieu et al. 
2007) and blinkers (Brooks ct al. 2004). Some recent observa- 
tions by EIS have also suggested that wc may be close to re- 
solving the cross field spatial scale of the fundamental coronal 
structures: see e.g. Warren et al. 2008a who derive loop filling 
factors of ^ 10%, or Tripathi et al. (2009) who obtained values 
larger than this except near the base of the loop they analyzed. 
This also sets an upper limit on the spatial scale of coronal heat- 
ing in loops, whether the fundamental structures are monolithic 
or composed of multiple threads. It does not mean, however, 
that the individual threads of the bundle are resolved. Based 
on magnetic fiux arguments, Priest et al. (2002) suggest that 
single Ti? ACE loops may be composed of 10 or more individual 
strands, and that the spatial scale we need to observe may be 
finer still. This is the case for nanoflare heating (Parker 1988) 
where the rcconncction takes place at the current shoots between 
threads (Klimchuk 2006). Current EUV and X-ray instrumen- 
tation are unable to resolve such structure, but the spatial reso- 
lution, excellent seeing, and stability of the Hinode Solar Opti- 
cal Telescope (SOT) allow it to observe at much higher spatial 
resolution than EIS or XRT. A direct comparison between chro- 
mospheric and coronal observations is of course difficult because 
of the large difference in instrumental spatial resolutions, and a 
lack of knowledge of the influence of the expansion of the field. 
Therefore, in this paper, we analyze only SOT observations of 
the magnetic field and chromospheric emission. We again study 
the core of the June 2007 region analyzed in Paper 1 and Paper 
2. Since the observations are new in themselves, we describe 
the magnetic field characteristics of the moss and active region 



core. Our main goal, however, is to study the variation of the 
magnetic flux, vector fleld, and chromospheric emission to try to 
uncover any behavior that could be related to the heating pro- 
cess and would allow us to set further constraints on the heating 
time-scale in this region. 

In §2 we describe the observations and data reduction proce- 
dures. In §3 we discuss the magnetic characteristics of the region 
and in §4 we introduce the modeling that we use to compare the 
time-scales of observational variability with typical loop cooling 
times. This is important because it is the crucial time-scale that 
differentiates between low and high frequency (effectively steady) 
impulsive heating. Wo thus put several statements about the 
evolution of the magnetic field on a quantitative basis. In §5 we 
examine the variability and evolution of the magnetic field and 
the chromospheric emission, and find that most of the activity 
takes place on time-scales longer than the computed theoretical 
loop cooling times. We also discuss departures from this picture. 
In §6 we discuss the spatial resolution issue once again and we 
derive the filling factor for the moss using EIS observations. The 
results suggest that SOT may be able to resolve the cross-field 
spatial scale of structures in the moss in individual pixels, so we 
further examine the magnetic flux and chromospheric variabil- 
ity on these size scales §7. The conclusions are presented and 
discussed in §8. 

2. Hinode AND TRACE Observations 

AR 10960 crossed the solar disk between 2007, May 30 and 
June 14. The region produced numerous C- and M- class flares 
during that period and was therefore the main observing target 
for most solar instrumentation on the ground and in space. In 
this paper we mostly use SOT data, but we also coalign and use 
TRACE 171 A filter images to identify the coronal features in 
the active region core, particularly the moss. An XRT Open/Ti- 
Poly image is also used for giving an overview of the hot core 
emission, and later we also use EIS data to determine the moss 
filling factor. 

The SOT is described in detail by Tsuneta et al. (2008). It 
consists of an optical telescope assembly (OTA, Suematsu et al. 
2008) that feeds a Filtergraph (FG) instrument and Spcctropo- 
larimeter (SP). The FG itself consists of a Broadband Filter 
Imager (BFI) and a tunable Lyot-type Narrowband Filter Im- 
ager (NFI). The NFI can obtain filtergrams, dopplergrams, and 
Stokes I, Q, U, and V images in a number of spectral lines formed 
ill the solar photosphere and chromosphere. The SP instrument 
obtains high-precision polarimetric scans in the Fe I 6301 A and 
6302 A spectral lines. The precision of the polarimetric calibra- 
tion is discussed in Ichimoto et al. (2008). As it is flown in space 
and has its own correlation tracker, SOT obtains high quality 
seeing-free and stable longitudinal or transverse magnetograms 
that can be used to diagnose magnetic field dynamics in the 
lower atmosphere. 

In this paper wc analyze both FG and SP data. Large FOV SP 
observations are used in §3 to describe the magnetic characteris- 
tics of AR 10960. These data were obtained from the SOT level- 
2 archive. As such they are outputs from the Milnc-Eddington 
gRid Linear Inversion Network (MERLIN) code developed at the 
Community Spectro-polarimetric Analysis Center (CSAC) at the 
High Altitude Observatory (HAO) by Bruce Lites and colleagues 
(Lites et al. 2007). MERLIN performs Levenberg-Marquardt 
least squares fitting of the full Stokes profiles obtained by the 
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Fig. 1.— Full Sun context images of AR 10960. Left Panel: SOHO/EIT 195 
magnetogram with the Hinode SOT FOV overlaid as a box. 

SP. Several assumptions about the atmospheric approximation 
and fit parameter initiahzation are made, for example, the source 
function is assumed to be linear with optical depth and the at- 
mosphere is assumed to be in local thermal equilibrium. The 
level-2 archive data are the parameters that give the best fit to 
the observed profile and much more detail is given on the level-2 
archive website. A key point of note is that no attempt is made 
to resolve the 180° azimuth ambiguity. This is discussed at the 
appropriate times below. 

The large FOV SP observations we analyze were obtained by 
scanning the slit over an area of 279.1" x 163.8" between 14:20:05 
and 15:23: 17UT on June 09. The resolution of this scan is 
0.30" X 0.32" per pixel and the time for individual polarimet- 
ric exposures was 1.6s. Exposures were obtained at each new 
scan position every 3.8s. In §5 we analyze a time-series of high 
cadence SP scans obtained between 21:56:03 and 22:59:31UT on 
June 08. The observing parameters were the same for this run, 
but the FOV was only 7"x512" so the time for each small scan 
is 28s. These data were also obtained from the level-2 archive. 

We also analyze a time-series of FG longitudinal magnetograms 
obtained in the Nal D 5896 A line at ~ 30s cadence. The spa- 
tial sampling of the Nal D Stokes images was 0.1 6" /pixels over a 
327.7" X 163.8" FOV with an effective exposure time of 14.2s cre- 
ated from polarimetric images with individual exposure times of 
0.12s. The time-series ran between 18:14:32 and 23:37:04UT on 
June 09, and the data were processed using the SolarSoft routine 
FG_PREP. The SOT observations were set to complement each 
other, with similar FOV scans with the SP and FG interspersed 
with high cadence time-series. They cover approximately 26 hr 
in total. 

To study the variability of chromospheric emission we also an- 
alyze a co-temporal FG time-series of Can 3896 A images taken 
between 18:14:54 and 23:37:26UT on June 09. Images were ob- 
tained at 30s cadence with an exposure time of 0.15s. The spa- 
tial sampling of these data is 0.109"/pixels over a 111.6" x 111.6" 



imago with the TRACE FOV overlaid as a box. Right Panel: SOHO/MDl 



FOV. To estimate the statistical noise in Can in §5.2 we use a 
very high cadence (8s) time-series of a plage region obtained on 
2007, February 20, between 11:57 and 13:37UT. All the Can 
data were processed using FG_PREP. 

TRACE also repeatedly observed AR 10960 during its passage 
across the disk, and movies of the region have been presented 
in Paper 1 and Paper 2. In §3 we use a 171 A filter image ob- 
tained at 14:37:23UT on June 09 to identify the moss areas. 
The FOV was 512" x 512" and the exposure time was 32.8s. We 
also use a 1600 A filter image obtained at 14:36:09UT to coalign 
with the SP data. The exposure time for this image was 0.86s. 
The data were processed and despiked using the SolarSoft rou- 
tine TRACE_PREP. In addition, we use an XRT Open/Ti-Poly 
image taken at 14:34:46UT with a 512" x 512" FOV and 0.06s 
exposure time. This image was processed using XRT_PREP. 

In section §6 we use EIS data to measure the moss filling fac- 
tor. The EIS instrument is described in detail by Culhane et al. 
(2007) and Korendyke et al. (2006). It observes with high spec- 
tral and spatial resolution (22.3mA and l" pixels, respectively) 
in short and long wavelength bands in the ranges 171-212 A and 
245-291 A. Several slits from 1" to 266" wide are available for 
making observations. In this paper we use a raster scan obtained 
with the observing sequence AR_velocity_map at 10:58:10UT on 
June 09. The sequence runs for about 5 hr 15 mins and takes a 
context image by stepping the 40" slit across a large FOV fol- 
lowed by a 1" slit scan over 330" x 304". We only use the 1" 
slit scan data here, the exposure time for which was 40s at each 
position. Many lines are included in the study, though we only 
need one density diagnostic line pair for our analysis. The data 
were processed using the default options in the SolarSoft routine 
EIS_PREP. Each of the lines used were then fitted with Gaussian 
profiles to obtain the line intensity. Further details are given in 
§6. 
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Fig. 2. — SP, TRACE, and XRT images of AR 10960. Top Panel: Hinode 
SOT/SP map corresponding to the FOV shown in Figure 1. The image 
shows the magnitude of Stokes V signal in the Fcl 6302 A profile. Middle 
Panel: Coahgncd TRACE 171 A filter image taken at 14:37UT during the 
SP scan and showing the moss emission. Approximately coaligned XRT 
image taken at 14:35UT during the SP scan and showing the location of 
the high temperature core emission. 





Fig. 3. — Coaligned SP and TRACE overlay, magnetic field inclination, 
and azimuth angle Top Panel: TRACE image overlaid in green on the 
SOT/SP map. The TRACE imago has been scaled logarithmically. The 
solid blue line outlines the major moss region and the dashed blue line is 
intended to exclude obvious loop emission. Middle Panel: f{0), where 9 is 
the inclination angle. Lower Panel: cos 20, where is the azimuth angle. 
The same contours from the upper panel arc overlaid in blue on the middle 
and lower panels. 



3. Magnetic Characteristics of AR 10960 

Here we describe the overall magnetic structure of AR 10960. 
To understand the relationship between the moss and the mag- 
netic field we need to coalign the TRACE image to the SP large 
FOV map. To do this, we first cut out the common area from the 
TRACE 1600 A image as determined from the FITS header co- 
ordinates. We then re-sampled the SP data to the lower TRACE 
resolution (note that all of the quantitative analysis is done on 
the original data), and then coaligned the 1600 A image with the 
SP measured non-directional magnetic field strength. From this 
procedure we corrected the uncertainty in the common area of 
the 1600 A image determined from the FITS header coordinates. 
We then extracted the correct area from the 1600 A image. 

This procedure worked well, but some stretching and rotation 
of the 1600 A image compared to the SP map was evident as 
a result of the difference in plate scale magnification and satel- 
lite orbital attitude. These effects were corrected after visual 
inspection. The rotation correction is approximately 1° counter- 
clockwise and the magnification correction is approximately 2% 
between the TRACE and SOT images. Furthermore, an E-W 
pixel shift of about 2 pixels was identified and removed. 

Having established the alignment between the SOT and 1600 A 
images, the coalignment with the 171 A TRACE image was fi- 
nally made by cross-correlating the two TRACE images after 



the 171 A image had been corrected for the inter- TRACE filter 
offsets reported by Handy et al. (1999). 

Figure 1 shows full disk images taken on June 09 by the 
SOHO Extreme ultraviolet Imaging Telescope (EIT, Delabou- 
diniere et al. 1995) and Michelson Doppler Imager (MDI, Scher- 
rer et al. 1995). The FOV of the TRACE data used to coalign 
with SOT is shown on the EIT image, and the FOV of the SP 
large FOV scan is shown on the MDI image. The extent of the 
active region and magnetic configuration are easily seen. Note 
the relative lack of complex emission or obscuration of the moss 
in the core of the region. This makes this particular region an 
ideal candidate for study. 

Figure 2 shows the magnitude of Stokes V signal measured 
in the Fel 6302 A line, a coaligned 171 A TRACE image, and 
an XRT Open/Ti-Poly image. The region is 13"/ by June 09, 
with well separated positive and negative polarity flux. There 
is strong moss emission above both polarities in the core, sepa- 
rated by a dark channel above weak magnetic flux. As pointed 
out previously by others (Katsukawa & Tsuneta 2005; Tripathi 
et al. 2008; Brooks et al. 2008; Title 2009), the plage and moss 
areas are predominantly unipolar with unmixed flux, though the 
spatial correlation between the fine detail is not clear (Berger 
et al. 1999a; de Pontieu et al. 2003). There is bright loop emis- 
sion to the south that obscures part of the moss. The XRT 
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Fig. 4. — Potential field extrapolation from the moss region indicated in 
Figure 2. Only 218 field lines are shown for clarity. 

image shows that hot loops surround this area with core emis- 
sion around the sunspots and covering the moss. The detailed 
relationship between warm and hot emission in this region was 
investigated in Paper 2. 

The results of the detailed coalignment between TRACE and 
the SOT/SP map are shown in Figure 3. An overlay of the 171 A 
image on the SP data is shown in the top panel. We identified 
the major moss emission visually with a contour level ~25% of 
the brightest loop emission in the 171 A image. This is shown by 
the solid blue line. Within this region, there is clearly bright loop 
emission. To exclude these structures, we drew another contour 
at the 45% level (dashed blue line). The pixel coordinates of 
all points inside the moss contour but outside the dashed region 
were recorded for all the images. From here on when we refer to 
the moss region, this is the area we are referring to. 

The magnetic field line inclination (9) and azimuth angle (</>) 
in the region are shown in the middle and lower panels of Figure 
3 with the same contours overlaid. In spherical coordinates, the 
field line inclination is along the line-of-sight so that an angle 
of 0° is directed towards the observer, and an angle of 180° is 
directed away from the observer. This will be relative to a radial 
field-line provided that the line-of-sight is normal to the surface, 
for example, when an active region is at disk center. For the 
data shown, the active region was at ~ 30° West of disk center. 

The inclination angle in the figure is represented by the function 

r e ife<^/2 
^(^)-\ Tv-e ife>^/2 W 

This image shows large values of the inclination as white, and 
smaller values as black. 

It has been pointed out previously that the magnetic field line 
orientation can be almost vertical to the solar surface (Kat- 
sukawa & Tsuneta 2005). The field rooted in the moss in this 
region appears to be only moderately inclined, with strongly in- 



clined field mostly around the edges of the moss or in the inter- 
moss lane. It is also notable that the inclination is coherent, i.e., 
there is no obvious mixing of widely differing inclination angles. 

The azimuth angle of the field is of course difficult to interpret 
because of the 180° ambiguity. For these data, the angle is mea- 
sured from to 180° from the Solar West position (RHS), but 
the value of the angle could be 180° in the opposite direction. 
Therefore, following Kubo et al. (2007), we show a cos 20 repre- 
sentation in Figure 3. This representation shows magnetic field 
oriented East- West as white, and field oriented North-South as 
black. 

The strong moss emission is characterized by field oriented 
East- West in both polarities. The orientation again seems coher- 
ent, with the field mostly changing direction around the edges 
of the moss. The edges then, would be the places where most of 
the shear in the magnetic field is located. The strongest shearing 
is close to the inter-moss region where the field is predominantly 
oriented North-South. 

4. Modeling 

It has previously been noted that the large scale pattern of 
TRACE 171 A brightness in the moss evolves slowly, with dy- 
namics of motion and variability mostly on small-scales (Berger 
et al. 1999a,b). Brooks et al. (2008) noted that although there 
was fine scale variability in the magnetic flux below the moss in 
the region they studied, the general pattern also evolved slowly. 
In this paper, we assess the stability and evolution of the mag- 
netic flux and vector field in AR 10960, but we also put our 
statements on a quantitative basis by comparing the observed 
time-scales to typical loop cooling time-scales computed from 
hydrodynamic simulations. Loop lengths vary throughout an ac- 
tive region, and since the radiative and conductive cooling times 
for a loop are dependent on the loop length, we first determined 
the distribution of lengths for loops relevant to the moss region 
using potential field extrapolations. The real magnetic field in 
the solar atmosphere is not likely to be current free (see e.g. De 
Rosa et al. 2009), however, the results from §3 show that the 
magnetic field is sheared mostly in the inter-moss region and 
that the field in the moss itself is unsheared with a small incli- 
nation angle. Since the moss is also unipolar, the field should 
escape almost vertically to the solar surface since it has nothing 
to connect to locally even if the field were non-potential. The 
potential field approximation, therefore, may be less problematic 
in the moss itself. In any case, the extrapolation is only used 
here to provide a realistic distribution of loop lengths. 

A magnetogram for extrapolation was prepared by weighting 
the SP magnetic field strength by the magnitude of the Stokes 
V signal. The weighting determines the polarity of the flux. An 
area of 208" x 160" around the core of the active region was then 
extracted. Field lines were computed for every pixel in the mag- 
netogram with a field strength between 5 and 5000G. Field lines 
rooted outside the moss region or that left the computational 
domain were then discarded. This resulted in a final dataset of 
23917 field lines. A subset of 218 are shown in Figure 4 overlaid 
on the magnetogram. 

An interesting result from this simulation is that the connec- 
tivity of the magnetic field in the moss is somewhat different 
than that of the quiet Sun. Close et al. (2003) showed that 
approximately 50% of the magnetic flux closes within 2.5Mm 
in the quiet Sun, with only 5-10% extending to heights greater 
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Fig. 5. — Rotated view of the potential field extrapolation of Figure 4 
showing thai, in contrast to the quiet Sun, most of the moss field connects 
to the corona. 

than 25Mm. As mentioned, the magnetic field in the moss has 

nothing to connect to locally in the unipolar regions, so it can 
extend higher. Figure 5 shows a side-view representation of our 
simulation with the height extension and flux closure percent- 
ages indicated. Considering only the field lines rooted in > 20G 
field (approximately comparable to the simulations of Close et al. 
(2003)), we find that about half extend beyond lOMm with a sig- 
nificant fraction (30%) extending to 25Mm and above. Many of 
the shortest field lines in our simulations are around the edges of 
the moss or crossing the inter-moss region. In these regions the 
potential field extrapolation is less likely to accurately represent 
the real sheared vector field, so the fraction of moss field lines 
that extend to significant heights is probably higher. This sug- 
gests that the moss is different than the quiet Sun in that most of 
the magnetized chromosphere is in fact connected to the corona. 
As we have shown (§3), this is also consistent with SP obser- 
vations that the departure from the radial field line direction is 
only moderate. 

The distribution of loop lengths for the full 23917 moss extrap- 
olated field lines is shown in Figure 6 (upper panel). The distri- 
bution shows significant numbers of short (< 30Mm) field lines 
and populations in the range 30-80Mm and 80-120Mm. The 
maximum length is 337Mm. To compute cooling times for the 
moss loops we prepared a grid of lengths spanning a wide range 
up to 260Mm for input into the hydrodynarnic code. For this 
simulation, we used the NRL Solar Flux Tube Model (SOLFTM) 
described in detail by Mariska (1987) and Mariska et al. (1989). 
Each loop on the grid was allowed to cool from a starting equi- 
librium temperature of ~ 4.8MK. The apex temperature at each 
computational time step is computed by averaging over the loop 
top. The loop is considered to have "cooled" when it reaches 
IMK. This is close to the formation temperatures of Fe ix and 
Fex, the spectral lines of which contribute significantly to the 
TRACE 171 A pass-band. The result of the simulation is shown 
in the lower panel of Figure 6. The cooling times range from 
300s for the shortest loops in the simulation, to 4800s for the 
longest loops. A linear fit to the results gives a relationship 

= -127.7 -M8.7L (2) 

where L is the loop length, and tc is the cooling time from 4.8 
to IMK. 
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Fig. 6. — Top Panel: Distribution of loop lengths from the potential field 
extrapolation for all closed field lines in the moss regions. Bottom Panel: 
Loop cooling time versus length calculated by SOLFTM for the distribu- 
tion of lengths in the top panel. Definitions are in the text. 



5. Variability and Evolution of Magnetic Field and 
Chromospheric Emission 

5.1. Magnetic Flux 

The time-scries of FG magnctograms taken on June 09 was 
used to assess the variability of the magnetic flux in AR 10960. 
For this global look, we co-registered the time-series using cross- 
correlation on rcsamplcd (half dimension) data. We then calcu- 
lated the average and standard deviation of the magnetic fiux 
in each pixel over the whole time-series; \B\ and a, respectively. 
We then calculated the quantity o-/\B\, that gives a measure of 
the variability in each pixel throughout the observations. This 
is over-plotted in rod on B in Figure 7 so that the locations of 
high and low variability in the region can easily be distinguished. 
Note that since B is averaged over the whole time-series, B high- 
lights the areas where the magnetic flux persists. Only the points 
were B is above the estimated statistical noise (see below) are 
plotted. 

As one would expect, the magnetic flux is persistent in and 
around the strongest field. It can also be seen that the variabil- 
ity is highest away from these regions. In particular, the vari- 
ability is low in the moss and high around the edges. It is also 
high around the neutral line that passes through the inter-moss 
lane. In previous studies it has been shown that most transient 
brightening activity in the EUV occurs around the neutral line 
of an active region, and that the magnetic flux pattern is persis- 
tent around the bases of high temperature loops (Brooks et al. 
2008). This plot is consistent with that picture. 

The movie (moviel.mpg) associated with Figure 7 shows the 
spatial distribution of the magnetic flux evolving slowly. To 
quantitatively asses this evolution we computed the linear Pear- 
son cross-correlation coefficients, r, between magnetograms sepa- 
rated by varying time-intervals. The Pearson coefficient is calcu- 
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Fig. 7. — Plot of (7/|i?| overlaid in rod on B for each pixel in the eoniplete time-series of SOT/FG data. The solid green line shows the area used for the 
cross-correlation results of Figure 8. Only the points where the magnetic flux is above the estimated statistical noise (16. 6G) are shown. A movie of the 
complete scries of magnetograms is available in the electronic edition of the manuscript as moviel.mpg. 



lated by dividing the covariance of two images by the product of 
their standard deviations, and was computed here for the boxed 
green area indicated in Figure 7 using the IDL routine CORRE- 
LATE. Since each magnetogram is separated by only 30s, r will 
be very high if the time-interval between them is small. As the 
separation is increased, the correlation will begin to break down 
because of the evolution of the spatial distribution of the flux. 
The data are interrupted by Hinode night time every orbit, so 
the cross-correlation is made only between magnetograms taken 
in the same orbit. 

Figure 8 shows the results. For magnetograms taken at 30s fre- 
quency r is close to 1.00 for the duration of the time-sequence, 
indicating a strong correlation between successive magnetograms 
as expected. As the separation between magnetograms is in- 
creased to 120s r falls, but is still maintained close to 0.98. With 
a separation of 960s r is close to 0.94, and with a separation of 
1920s r is maintained close to 0.90. 1920s is equivalent to the 
cooling time for a loop of length llOMm in our hydrodynamic 
simulations. Since 85% of the loops in our fleld extrapolation are 
shorter than this they cool from ~ 5 to 1 MK on shorter time- 
scales. This indicates that the spatial distribution of the mag- 
netic flux maintains a strong correlation for time-scales longer 
than the cooling times for at least 85% of the loops extrapolated 
from the moss regions. With the data at hand, it is not possible 
to assess the pattern of magnetic flux over longer time-scales. 

The variation of the magnetic flux in small areas as a func- 
tion of time was then investigated. The analysis is not sophis- 
ticated. We used the full resolution data and extracted a 1000 
pixel'^ area from each magnetogram. This time-series was then 
co-registered via cross-correlation of successive magnetograms. 
Figure 9 shows one of the FG magnetograms in the time-series. 
Eleven small boxes are overplotted scattered throughout the pos- 
itive and negative polarity moss regions. One in the inter-moss 



region is also included. The boxes are 2"x2". This size was 
chosen because it is comparable to the spatial resolution of EIS, 
which was measured to be close to 2" in the laboratory pre- 
launch (Korendyke et al. 2006). It is about the same size as a 
low resolution MDI pixel. 

In order to understand whether the variations in small areas are 
signiflcant or not, we need to estimate the statistical uncertainty. 
Since the FG does not do a full magnetic vector inversion we can- 
not quantify the various sources of noise. Therefore, we adopted 
the following approximate method. We selected a relatively quiet 
region, shown by the large box in Fig 9, and formed histograms 
of the difference in magnetic flux per unit area between succes- 
sive magnetograms in the time-series. These histograms were 
then fitted with Gaussian functions and the standard deviation 
(cr) measured. The standard deviations are plotted as a func- 
tion of time in the upper left panel of Fig 10. The average value 
is 16. 6G with ~ 15% variation in the standard deviations as a 
function of time. We adopt this average value as our estimate of 
the statistical noise. There may be signals below this value, but 
they would be difficult to reliably distinguish from noise, with 
the caveat that the quiet region selected is in an active region 
so is only 'relatively' quiet. It is likely that the true statistical 
noise level is lower. This estimate should be taken into account 
when considering the results for the moss and inter-moss regions 
below. 

Figure 10 shows the results of this analysis. The interpreta- 
tion is complicated because it is difficult to distinguish changes 
in the magnetic flux due to features evolving, or moving in or 
out of the field of view. As noted by Brooks et al. (2008), the 
constancy of the spatial distribution of magnetic flux in an ac- 
tive region can be preserved by magnetic features moving along 
similar paths, i.e., the field is dynamic, but the pattern is main- 
tained. We are interested, of course, in the evolution of features, 
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Fig. 8. — Variation of the linear Pearson correlation coefficient between magnetograms separated by varying time-intervals as a function of time. The 
correlation is made between time-series with magnetograms taken at frequencies of 30s, 120s, 960s, and 1920s. The last value is equivalent to the cooling 
time for a loop of length llOMm. 



but also want to eliminate cases where wo end up tracking some- 
thing else. We have therefore made an effort to select regions 
where obvious motions in or out of the box are reduced. The 
movie (movie2.mpg) associated with Figure 9 allows the reader 
to independently judge this analysis. 

The positive polarity moss box in the top row of Figure 10 
shows am average magnetic field of 510G and around 22% vari- 
ation over the duration of the observations (> 5 hr). This is 
clearly much longer than the theoretical cooling times. The neg- 
ative polarity moss box in the top row shows a comparable av- 
erage magnetic field of -490G and ~ 27% variation during the 
observations. The majority of the boxes in Figure 10 show fluc- 
tuations around this level of 15-30%. These values for the moss 
arc a little higher than the results found for coronal intensities 
and velocities in papers 1 and 2. A few cases show larger varia- 
tions at the 40-50% level, but it is clear that this also reflects the 
fact that the magnetic flux is evolving slowly. Note, for example, 
the box in the lower right panel of Figure 10. This box shows a 
variation of ~45% during the observing period. This is a result, 
however, of a slow evolution from around -550G at 18:30UT to 
around -200G at 23UT rather than a sudden change. In fact, 
even in this box, a less than 30% variability is maintained if we 
consider just the first, or last, 3 hours of observations. If we con- 
sider only the first hour of the observations the magnetic flux in 
all of the boxes varies by less than 10%. This is still considerably 
longer than most of the simulated loop cooling times: only loops 
longer than ~200Mm persist for an hour, and only about 2% of 
the loops in our simulation are this long. 

The inter-moss region shows an average magnetic fleld of 15G 
and the highest variation of all the boxes (~ 60%) This is con- 



sistent with our expectations since the inter-moss region was al- 
ready identified as being more variable in Figure 7. As pointed 
out earlier, transient EUV brightenings have been found to be 
preferentially located around the neutral lines in other active 
regions. The inter-moss region also seems to be where the ma- 
jority of the flares and transients occur in this region (see, e.g., 
the TRACE movies in papers 1 and 2). 

Some areas of the moss do show evidence of dynamic changes 
in magnetic flux similar to that in the inter-moss box. Also, it 
is unclear how significant 15-30% variations in moss magnetic 
flux are over time-scales of several hours. Variations comparable 
in magnitude do occur on shorter time-scales (tens of mins). 
These changes are on the order of 100-200G and are statistically 
significant. Therefore, they could be taken as evidence for low- 
frequency impulsive heating. Note that changes on 20-30 mins 
time-scales correspond to the cooling times for loops in the 70- 
lOOMm range. We find, however, that only a small fraction 
(~10%) of our simulated moss loops would be expected to cool 
on these time-scales. This suggests that such heating is not 
significant here, and may also be consistent with the lack of 
"warm" EUV loops cooling in the core of this region. 

The 100-200G magnetic fiux changes correspond to ~ 10^**Mx 
in the boxed regions, and would be equivalent to smaller than 
micro-flare size energy releases in the corona. Further work is 
needed to see if these could be related to, or signatures of, the 
flare-like intensity variations intermittently seen in hot X-ray 
core loops (Warren ct al. 2007) or hot transient events (Shimizu 
1995). Here we focus on determining whether such variations 
show evidence of propogating energy into the chromosphere. 
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5.2. Can H Intensity 

The variation of chromospheric emission in the Can 3896 A 
hne as a function of time was investigated using the time-series 
of images taken on June 09. The time-series was co-registered 
to the first image by cross-correlating successive images. Figure 
11 shows an example image taken at 18:14:55UT with the same 
eleven small boxes and large control box as in Figure 9 overlaid. 
The placement of these boxes was made by determining the inter- 
FG offsets and scalings by co-registering this Can image with a 
near-simultaneous Nal D magnetogram. This was achieved by 
cross-correlating a common area from the two images and re- 
sampling the magnetogram to the higher resolution of the Call 
image. 

As with the magnetic flux analysis, we need to estimate the 
statistical uncertainty in order to gauge whether variations in 
the Can emission are significant. For this purpose we analysed 
a time-series of very high cadence (8s) Call images of a plage 
region obtained between 11:57 and 13:37UT on 2007, February 
20, with an exposure time of 0.41s. The data were co-registered 
by cross-correlating successive images and intensity difference 
images for the entire time-series were prepared. Histograms of 
the intensity differences over the full FOV (56" x 28") were then 
formed for each image in the time-series and fit with Gaussian 
functions to determine the standard deviation (a) . The standard 
deviations are plotted as a function of time in the upper left panel 
of Fig 12. The average value is 9DN/s with only 4% variation 
in the standard deviations as a function of time. We adopt this 
average value as our estimate of the statistical noise. 

Figure 12 also shows the variation of intensity in the moss 
and inter-moss boxes. An animation (movie3.mpg) shows the 
locations of the boxes, and their stability, during the time-series. 
Remarkably, the average intensities in all the boxes (including 
the inter-moss box) are within 10% of each other. Furthermore, 
the intensity variation is less than 10% in all the boxes for the 
duration of the observations (> 5 hours). The largest variation 
can be seen in the lower left hand box. As with the magnetic 
flux, this is again clearly a slow evolution from ~ 1200DN/s at 



the start of the sequence to ~ 900DN/s at the end. Comparable 
magnitude variations of 100-200DN/s are again seen on time- 
scales of tens of minutes, and these are significantly above the 
noise estimate. 

We computed the correlation coefficients between the Call in- 
tensities and magnetic fluxes for each of the boxes. Intriguingly, 
about half show a correlation coefflcient > 0.5, indicating a weak 
positive correlation and suggestive that the 100-200G changes 
in the magnetic flux do propogate energy and heating that leads 
to emission that is detectable above the noise level in the chro- 
mosphere, though it is not universally seen. Sakamoto et al. 
(2008) have recently also reported the detection of fluctuations 
in TRACE 171 A observations above the photon noise level that 
have durations that agree well with their estimated loop cool- 
ing times, albeit for a different active region. Our analysis of 
EIS 195.119 A observations of this region (paper 1) showed only 
~ 15% variations in intensity over many hours, but a detailed 
comparison with the noise was not made. The lack of a con- 
sistent detectable signature in the chromosphere and a possible 
signature in the TRA CE data could indicate that the energy pro- 
pogation is amplified in the corona, or is in fact released there. 
Further work on this issue is warranted, including comparisons of 
intensity variations with numerical models ( Antolin et al. 2008) . 

5.3. Magnetic Field Vector 

Figure 13 shows the same magnetic field inclination map rep- 
resented by f{9) as in Figure 3 and defined in Equation 1. The 
FOV of the narrow (7"x512") SP slit scan from June 08 is over- 
laid. The FOV has been rotated to the time of the large FOV 
scan (14:20 on June 09). The same 171 A contours as before are 
shown, ft can be seen that the narrow scan crosses the inter-moss 
region South-North and reaches both positive and negative po- 
larity moss regions. Both regions also appear to be moderately 
inclined to the line-of-sight (black) whereas inter-moss region 
shows strong inclination (white). Two regions are selected in the 
moss (the small boxes within the FOV in Figure 13) to study 
the variation of the field line inclination and azimuth angle. The 
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Fig. 10. — Variation of moss magnetic field - flux per unit area - in Gauss (G) as a function of time. Top left: standard deviation of liistograms of the 
pixel-to-pixel difference between successive magnetograms in the control box shown in Figure 9. Other panels: evolution of the magnetic field averaged over 
the small boxes in Figure 9. The crosses show the individual points, and the thick solid lines show detrended 5 mins averaged curves to de-emphasize the 
scatter. 



areas of the boxes are again 2" x 2" . 




Fig. 11. — FG Can 3896A magnotogram taken at 18:14:54UT on 2007, 
June 09. Eleven small boxes are overlaid and the variability of the Call 
intensity in these regions as a function of time is plotted in Figure 12. The 
moss contours from Figure 3 are overplotted in blue. 

Figure 14 shows the results. Consistent with the measurements 
in §3 both moss regions show moderate inchnation from the hne 
of sight with an average value of 27° . The variation is less than 
15% over the time-series (21:56:03-22:59:31UT). The azimuth 



angle changes from ~ 120° at the start of the time-series to ~ 90° 
by the end, but as with other quantities, this is a slow evolution, 
and the variation from the average (~ 109°) is less than 15% 
over the whole period. Note that this time-series lasts for 3808s, 
which is again much longer than the theoretical cooling times for 
loops extrapolated from the moss regions calculated in §4. 

One criticism of such simple analyses is that variability on small 
scales is smoothed out by the averaging over large areas. It is 
possible that coronal emission is quasi-steady over many hours on 
size scales at the limit of current EUV/X-ray instrumentation, 
but that future instrumentation will find higher variability on 
smaller spatial scales. What tends to happen is that quantities 
appear relatively steady for many hours, within some threshold, 
when averaged over large areas. These quantities also appear 
steady when averaged over smaller areas, but for shorter periods 
of time, though the duration can be extended if the threshold is 
relaxed. For example, the most variable quantity in Figure 14 
is the inclination angle in the positive polarity box. It varies by 
less than 15% over the observations period. Individual pixels in 
this box, however, have larger variabilities: ~ 34% on average. 
Relaxing our definitions slightly, however, we are still able to say 
that > 80% of the individual pixels vary by less than 30% over 
time-scales comparable to the cooling times for 85% of the loops 
in our simulations (1920s). 

The key question of course is how these time-scales, spatial 
scales, and thresholds, compare to the critical values for coronal 
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Fig. 12. — Variation of moss Call intensity (DN/s) as a function of time. Top loft: standard deviation of histograms of the pixcl-to-pixel difference between 
successive images in the control box shown in Figure 11. Other panels: evolution of the Call inensity averaged over the small boxes in Figure 11. The 
crosses show the individual points, and the thick solid lines show detrended 5 mins averaged curves to de-emphasize the scatter. 



heating. We have shown that the chromospheric Can emission 
and all the magnetic activity in AR 10960 evolves slowly, and 
our quantitative analysis also indicates that the variability in 
magnitude is only at a low level on time-scales longer than the- 
oretical loop cooling times in small 2"x2" boxes. The movies 
presented in this paper clearly suggest, however, that there is a 
lot of small-scale dynamic activity. The time-scale for loop cool- 
ing is clearly resolved in these observations, but the important 
spatial scales may not be. Next we investigate the moss filling 
factor in this region using EIS observations in order to at least 
indirectly infer the size scale of a fundamental loop envelope. 

6. EIS Filling Factor Measurements 

Previous studies have shown that intensities in the moss scale 
linearly with the loop base pressure and independently of the 
loop length (Martens et al. 2000; Vourlidas et al. 2001). As 
pointed out by Warren et al. (2008b), this means that any phys- 
ical model that yields the same base loop pressure will predict 
the same moss intensity regardless of the loop length. Warren 
et al. (2008b) calculated full solutions to the hydrodynamic equa- 
tions, assuming steady heating, for a grid of loop lengths from 
10-lOOMm. This grid sufficiently covers most of the distribution 
of loop lengths we computed for the moss in §4. For each loop 
solution, the density and temperature around the loop is known 
and the intensity at each position for any EIS line can be cal- 
culated. By integrating the intensity over the lower 5Mm of the 
loop the footpoint intensity can be computed. The intensities for 



the grid of loops scale linearly with pressure and can be fitted 
with a function of the form Ix = clPq, where / is the intensity 
and Po is the base pressure. Warren et al. (2008b) state that the 
linear relationship between intensity and pressure breaks down 
at low pressures, so their fits are restricted to values above log 
P, = 16 cm^^* K. 

The intensity ratios of selected lines can also be directly re- 
lated to the simulated base pressure. By determining the moss 
pressure from observed intensity ratios, a simulated intensity for 
an individual line can be calculated from the fit and compared 
to the observed intensity for that line. The filling factor is then 
introduced (if needed) to bring the simulated and observed in- 
tensities into agreement. We followed this procedure here for 
AR 10960. 

Following the analysis in paper 2 we use a density threshold to 
identify the moss pixels. Several density diagnostics are avail- 
able in the EIS wavelength bands and a detailed discussion of 
comparisons between them has been presented in Young et al. 
(2009). In paper 2 we discussed these comparisons and decided 
to use the Fexill 202.044/203.826 A ratio because of the con- 
sistency between densities derived from this ratio and densities 
derived from Six ratios in small active regions, bright points, 
and the quiet Sun i.e. regions where the density is lower than in 
the moss and the Six lines are sensitive. For comparison with 
that work, we use this ratio again here. The analysis is applied 
to the EIS slit raster scan taken at 10:58:10UT on June 09, and 
discussed in §2. The FeXlll 202.044 A line was fitted at every 
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Fig. 14. — Variation of magnetic field inclination and azimuth angle as a 
function of time for the positive and negative polarity moss regions shown 
as boxes in Figure 13. The inclination is represented by f(0) (see text). 
Less than 15% variation is seen in all quantities over a period of about 1 
hr. 

blend of two lines and is further blended at 203.734 A with an 
Fexil line (Brown et al. 2008). We performed a double Gaussian 
fit to this feature at every pixel to remove the Fexil blend. To 
do this, we fixed the separation of the two components at 0.1 A 
and forced them to have the same width. 

Images of the active region core formed from the line fits are 
shown in Figure 15. The full raster spanned several passes 
through Hinode night, so the images only show a central area 



of 80" X 140". Electron densities were calculated for these re- 
gions using the measured line ratios and the CHIANTI database 
version 6.0 (Dere et al. 1997, 2009). A contour set at 40% of 
the maximum intensity is drawn on the Fexill 203.826 A im- 
age of Figure 15 to highlight the moss regions. Pixels within 
these contours were selected and identified as moss if the calcu- 
lated electron density exceeded log Ne = 9.6 cm~^. These are 
shown as crosses on the Fexiii 203.826 A image. For each of the 
identified pixels, the observed intensity ratio was computed after 
subtracting contaminant background emission. This background 
intensity was measured in the small box in the inter-moss region 
on the Fexill 203.826 A image. The new ratios were then used 
to derive the base pressure by interpolation from the simulated 
grid of solutions. Figure 16 shows the values for the moss pixels 
overplotted on the theoretical line ratio vs base pressure curve. 
Although several moss pixels approach the high pressure limit 
of the ratio, it can be seen that the vast majority of them fall 
in the sensitive range of the curve. The minimum base pressure 
in the moss is log Pe = 16.3 cm~^ K, which is higher than the 
lower limit of the power law fits of Warren et al. (2008b). With 
the base pressure established, the Fexill 203.826 A intensity was 
simulated using the fit for Fexill 203.826 A. The coefficients of 
the fit in this case are a — -11.94 and b = 0.99. 

As expected, the intensities thus simulated are much higher 
than observed, so a filling factor needs to be introduced to bring 
them into agreement. The distribution of filling factors for the 
moss pixels is shown in Figure 17. The majority of values fall in 
the 10-20% range with the median value being ~16%. This is 
in agreement with previous measurements of moss filling factors 
(Fletcher & de Pontieu 1999; Warren et al. 2008b). 

These results suggest that the fundamental size scale of struc- 
tures in the moss are about 10-20% of the EIS spatial resolution. 
This is consistent with similar measurements of coronal loops 
with EIS based purely on observational factors (Warren et al. 
2008a), and also sets an upper limit on the cross-field scale of 
coronal heating, which will be much smaller if loops are com- 
posed of multiple threads. Note that the filling factor is an area 
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Fig. 15. — EIS images of the core of AR 10960 formed from fits to tire den- 
sity sensitive Fexill lines. The moss regions are highlighted by a contour 
at 40% maximum intensity. The small white box shows the inter-moss 
region used for background subtraction, and the crosses show the moss 
pixels with log Ng > 9.6 cm^^ used for determining the filling factor. 



filling factor. EIS has 1" spatial pixels, so taken at face value a 
next generation instrument needs 0.3-0.4" spatial pixels to re- 
solve the median moss filling factor. Interestingly, this is larger 
than the size of the SOT FG and SP pixel scales. Although EIS 
(and XRT) cannot directly resolve these scales, SOT can resolve 
these scales below the moss. 

7. Magnetic flux evolution and chromospheric 
emission in individual pixels 

Having inferred the apparent cross-field size scale of funda- 
mental structures in the moss and established that the FG can 
resolve this scale in individual pixels, we revisit our analysis of 
the variability in magnetic flux and chromospheric emission in 
the moss region (§5). Previously we showed that the average 
magnetic field in most of the boxes varied by ~ 15-30% while 
the chromospheric Can 3896 A intensity varied by ~10%. Note 
that these changes include changes due to motions in and out 
of the field of view as well as variations in magnitude. These 
motions may have a more pronounced effect on the results for 
individual pixels because the features can traverse the pixel size 
scale faster. It is worth pointing out that features moving in 
and out of a pixel can have the effect of canceling each other to 
some degree so that variability may be reduced. This would also 
be true if the fundamental structures are composed of multiple 
threads. 

Figure 18 shows the distribution of percentage variabilities (ra- 
tio of standard deviation a to average magnetic field B) for all 
the pixels in two of the boxed regions. Each box contains 156 
pixels. The values are calculated for a time-period of 1920s. 
Recall that 85% of the loops in our simulation cool on shorter 
time-scales than this. In both cases > 75% of the pixels show 
variabilities less than 20% and the median values for both dis- 
tributions are ~ 13%. The variability is clearly greater in indi- 
vidual pixels over the full observations period, but this is also 
partly attributable to the lifetime of individual features. For the 
record, more than 2/3 of the pixels in the boxes show less than 
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Fig. 16. — Theoretical line ratio versus base pressure curve for the grid of 
loop models discussed in the text. The values for the moss pixels high- 
lighted in Figure 15 are overplotted as crosses. 



45% variability for the whole > 5 hr period. 

Figure 19 shows the distribution of percentage variabilities of 
the Call 3896 A intensities for all the pixels in the same two 
boxed regions, computed for the same time-period. In this case, 
each box contains 306 pixels. In both cases ~ 85% of the pixels 
show variabilities of less than 10% and the median values for 
both distributions are ~ 6%. Again, the variability is larger 
over the full observations period but, at least for the selected 
boxes, all the pixels show less than 15% variability for the whole 
> 5 hr period. 

To check whether the two examples are representative, we show 
in Figure 20 the distribution of percentage variabilities of the 
Ca II 3896 A intensities and unsigned magnetic flux for all indi- 
vidual pixels within all the moss boxes. The Can distribution 
has a median around 6% with 99% of the pixels showing variabil- 
ities of less than 15%. The median of the magnetic flux distri- 
bution is around 15% with more than 2/3 of the pixels showing 
variabilities of less than 20%. These values were again computed 
for a time-period of 1920s. As discussed in §5, the magnetic flux 
in the moss shows greater variation, but this does not seem to 
register strongly in the chromosphere. 

8. Discussion 

We have described the characteristics and evolution of the 
longitudinal and vector magnetic field and chromospheric Call 
emission in the core of an active region observed by Hmode and 
TRACE. Consistent with previous studies, we found that the 
moss in this region is unipolar, the field lines are only moderately 
inclined, and the spatial distribution of magnetic flux in the core 
evolves slowly. We put the last statement on a quantitative basis 
by comparing the evolution time-scales with the theoretical cool- 
ing times computed from hydrodynamic simulations of coronal 
loops extrapolated from SP magnetograms of the moss regions. 
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Fig. 17. — Distribution of filling factors computed for the moss pixels. 



These extrapolations also showed that the magnetic connectivity 
in the moss is different than in the quiet Sun because most of 
the magnetic field extends to significant coronal heights (there 
is nothing to connect to locally in the unipolar regions). We 
also showed that the field line inclination and horizontal compo- 
nent are coherent in the moss, and most of the shearing of the 
field occurs in the inter-moss region, or at the edges of the moss, 
where the magnetic flux variability was shown to be greatest. 
The magnetic vector, flux, and Can intensity in 2"x2" boxes in 
the moss also do not show signiflcant variability on time-scales 
longer than loop cooling time-scales when averaged over these 
areas. Though flux changes on the order of 100-200G and Ca ii 
intensity changes on the order of 100-200DN/s are observed on 
20- 30 mins time-scales, we find only weak evidence that these 
flux changes consistently register in the chromosphere. We also 
determined the filling factor in the moss from EIS observations: 
~16%. This is consistent with previous studies, and suggestive 
that the cross-field scale of the fundamental structures in the 
moss are larger than the size of an individual SOT FG pixel 
(0.16"). The mean variabilities of the magnetic flux and chro- 
mospheric emission in individual FG pixels was also found to 
be ~15% and ~6%, respectively, on time-scales longer than the 
computed loop cooling times. 

EUV and X-ray observations of the moss in the core of this 
region have already suggested that the heating is steady or ef- 
fectively steady i.e. heating events occur with a rapid repetition 
rate. The results presented here show that the magnetic fleld 
and chromospheric emission also evolve slowly, and remain rela- 
tively steady even at the highest spatial resolution we have ever 
observed. They could therefore be interpreted as further evi- 
dence supporting the quasi-steady picture. The short time-scale 
100-200G variations we see could be interpreted as evidence of 
low frequency impulsive heating. We flnd, however, that only a 
small fraction of our simulated loops (10%) would be expected 
to cool on these time-scales. 



Fig. 18. — Distributions of magnetic field variability (a/B) expressed as 
percentages for the individual pixels within a positive and negative polarity 
moss box from Figure 9. 



If the heating events do occur frequently and are rcconnection 
related, one would also expect to see changes in the magnetic 
field on a comparable time scale. Further studies of the relation- 
ship between the detailed small scale changes in the magnetic 
flux that we see and variations in the EUV/X-ray intensities are 
needed to address this issue. At present it is unclear theoreti- 
cally how important changes of 10-20% are in terms of energy 
input. The magnetic field inclination, for example, is indica- 
tive of the horizontal stress component of the field and could 
be directly realted to the heating. Small variations of this an- 
gle do not necessarily imply that the heating variations are also 
small. One nuance is that the measurements of changes in the 
magnetic field/flux are a combination of changes due to trans- 
lational motions as well as variations in magnitude, so direct 
conversion into energy is not unambiguous. Nevertheless, the 
quiet Sun chromospheric and coronal heating requirements are 
~4xl0® erg cm~^ s~^ and ~3xl0^ erg cm~^ s"'^, respectively 
(Withbroe & Noyes 1977; Aschwanden 2004), suggesting that 
10% changes in chromospheric output could be enough to sus- 
tain the quiet corona. The situation is more complex, however, 
in active region loops. Kano & Tsuneta (1996) and Katsukawa 
& Tsuneta (2005) quoted estimates of ~10^ erg cm~^ s~^ for 
'warm' loops and ^lO'^ erg cm~^ s~^ for hot loops. This com- 



erg cm s~^ 

3^ erg cm~^ s~^ 
pares to ~2xl0^ erg cm~^ s~^ for the chromosphere in active 
regions (Withbroe & Noyes 1977; Aschwanden 2004). Assuming 
that the strong basal levels of magnetic flux and Ca ii emission in 
AR 10960 match this heating requirement, these numbers sug- 
gest that the 10-20% changes we observe on 'warm' EUV loop 
cooling time-scales of 20 30 mins could be enough to heat these 
loops impulsively at low frequency. These changes, and the 10- 
20% variations that occur over time-scales that are longer than 
a loop cooling time, would not be enough, however, to heat the 
hot loops. 

Katsukawa & Tsuneta (2005) estimated the energy flux into 
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Fig. 19. — Distributions of Can intensity variability (cr/T) expressed as 
percentages for the individual pixels within the same two moss boxes shown 
in Figure 18. 



Fig. 20. — Distributions of Can intensity variability (t/T) and magnetic 
field variability {a/B) expressed as percentages for all the pixels in all the 
moss boxes shown in Figure 9. 



coronal loops as a result of the dissipation of magnetic energy 
built up in coronal current sheets due to the braiding of the field 
by photospheric motions. They found that for a moss magnetic 
field strength measurement of 1.2kG, the heating requirement 
of hot loops (10^ erg cm^^ s^^) could be met. In their esti- 
mates, however, the energy flux is related to the magnetic flux 
as F oc i?^, so that once again a 10-20% change in magnetic 
flux may meet the heating requirements of the 'warm' loops, but 
docs not appear to supply enough energy to heat the hot loops. 
A definitive statement on this issue awaits more detailed theo- 
retical modeling. These results suggest, however, that only the 
continuous input from the strong basal levels of magnetic flux 
and chromospheric emission, not the variations, appear to be 
enough to power the hot core loops rooted in the moss. 

The results presented here, together with those in Paper 1 and 
Paper 2 and Antiochos et al. (2003), show a lack of variability 
in many types of diagnostic signatures: magnetic flux observa- 
tions by SOT, coronal velocity measurements by EIS, EUV and 
X-ray intensities from TRACE and XRT, etc., and constitute a 
compelling body of evidence. Future instrumentation, however, 
may reveal all of these quasi-steady properties to be a distrac- 
tion. It is possible, for example, that the magnetic field could 
appear steady in the lower atmosphere while it is braided in the 
corona and energy is released there. Heating by impulsive events 
could occur at coronal heights with no visible signature in the 
magnetic field at lower heights. The lack of variability in the 
EUV and X-ray data could indicate that the impulsive events 
occur on sub-resolution scales. Low frequency impulsive heating 
on unresolved scales could give the impression of high frequency 
(effectively steady) heating. These suggestions, of course, are 
difficult to test, since we do not know the spatial scale of the 
heating. Priest et al. (2002) suggest that the size scale of the 
fundamental kG fiux tubes in the photosphere have diameters 
of 100km, and based on these arguments Klimchuk (2006) sug- 



gests that a comparable spatial resolution may be able to resolve 
the strands making up the corresponding loops, though the cur- 
rent sheet interface may be smaller. We have shown here that 
the magnetic flux and chromospheric emission are quasi-steady 
on spatial scales approaching 100km. This is 4000 times smaller 
than studied by Antiochos et al. (2003) and should rule out mod- 
els that predict low frequency impulsive heating on larger scales. 
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